Microscopic self-energy calculations and dispersive optical-model potentials 
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Nucleon self-energies for 40,48,60 Ca isotopes are generated with the microscopic Faddeev-random- 
phase approximation (FRPA). These self-energies are compared with potentials from the dispersive 
optical model (DOM) that were obtained from fitting elastic-scattering and bound-state data for 
40,48 Ca. The ab initio FRPA is capable of explaining many features of the empirical DOM poten- 
tials including their nucleon asymmetry dependence. The comparison furthermore provides several 
suggestions to improve the functional form of the DOM potentials, including among others the ex- 
ploration of parity and angular momentum dependence. The non-locality of the FRPA imaginary 
self-energy, illustrated by a substantial orbital angular momentum dependence, suggests that future 
DOM fits should consider this feature explicitly. The roles of the nucleon-nucleon tensor force and 
charge-exchange component in generating the asymmetry dependence of the FPRA self-energies are 
explored. The global features of the FRPA self-energies are not strongly dependent on the choice of 
realistic nucleon-nucleon interaction. 
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I. INTRODUCTION 

The properties of nucleons propagating in the nucleus 
exhibit characteristic deviations from the naive shell- 
model picture. At positive energies, corresponding to the 
domain of elastic scattering, there is unambiguous evi- 
dence that the potential (or self-energy) that a nucleon 
experiences is absorptive [lH3|. This simple observation 
has important implications, since it shows that nuclear 
states cannot be interpreted only in terms of a simple 
shell-model potential that is real and independent of en- 
ergy. The importance of the dynamic aspects of the nu- 
clear shell model was recognized in Ref. j4| ■ The link be- 
tween optical potentials and the traditional bound-state 
shell model was explored by Mahaux and Sartor H and 
extensively reviewed in Ref. @. These authors realized 
that information on nucleon propagation at positive en- 
ergy influences the properties of the real nuclear potential 
at negative energy, since the nucleon self-energy obeys a 
dispersion relation that links the real part to its imagi- 
nary part at all energies (see e.g. Ref. Q)- Mahaux and 
Sartor exploited standard representations of the imagi- 
nary part of the optical potential in terms of volume and 
surface contributions. They further assumed that the be- 
havior of the imaginary potential was similar near both 
sides of the Fermi energy and used a subtracted form of 
the dispersion relation to obtain the corresponding real 
part. By performing this subtraction at the Fermi en- 
ergy, only the additional knowledge of the real potential 
at that energy is required. The resulting optical potential 
is now called the dispersive optical model (DOM). 

Recent applications of the DOM have concentrated on 
the nucleon asymmetry dependence by simultaneously 
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fitting data pertaining to different calcium isotopes [8|, |9| 
and to spherical isotopes up to Tin and 208 Pb Such 
an analysis can be utilized to predict properties of iso- 
topes with larger nucleon asymmetry by extrapolating 
the DOM potentials. Such data-driven extrapolations 
present a reliable strategy to approach and predict prop- 
erties of isotopes towards the respective drip lines, since 
they can be tested by performing corresponding experi- 
ments. An important feature extracted from this anal- 
ysis is the increase in surface absorption of protons for 
increasing nucleon asymmetry. While this trend is unam- 
biguous, there is no clear understanding of the underlying 
dynamics responsible for it. A much weaker and opposite 
trend was inferred for neutrons. It is therefore useful to 
study the nucleon self-energy — which is the microscopic 
counterpart of the DOM — to clarify this behavior and 
provide a deeper understanding of the DOM potentials. 

The Green's function method is ideally suited to pur- 
sue a microscopic understanding of the nucleon self- 
energy at both positive and negative energies [ll|. The 
most sophisticated implementation of the Green's func- 
tion method considers the role of long-range or low- 
energy correlations in which nucleons couple to low-lying 
collective states and giant resonances. This is accom- 
plished by using the random phase approximation (RPA) 
to calculate phonons of particle-particle (hole-hole) and 
particle-hole type. These are then summed to all orders 
in a Faddccv summation for both two-particle-one-hole 
(2plh) and two-hole-one-particle (2hlp) propagation. 
This approach is referred to as Faddeev random phase 
approximation (FRPA) [12, [H| ■ This method is size ex- 
tensive and has been successfully benchmarked for soft 
interactions in purely ab-initio calculations for 4 He (3], 
giving results of comparable accuracy to coupled-cluster 
theory 

The FRPA was originally developed to describe the 
self-energy of the double closed-shell nucleus 16 O [HI, 
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Il5j . The method has later been applied to atoms and 
molecules (HOS and recently to 56 Ni [Hi and 48 Ca [3. 
The ab initio results of Ref. [l8[ are in good agreement 
with (e, e'p) data for spectroscopic factors from Ref. [lj| 
and also show that the configuration space needed for the 
incorporation of long-range (surface) correlations is much 
larger than the space that can be utilized in large-scale 
shell- model diagonalizations. In Ref. [20[, the FRPA was 
employed to calculate proton scattering on 16 O and ob- 
tain results for phase shifts and low-lying states in 17 F. 
However, the properties of the self-energy at larger scat- 
tering energies which are now of great interest for the 
developments of DOM potentials was not addressed. In 
particular, one may expect to extract useful information 
regarding the functional form of the DOM from a study 
of the self-energy for a sequence of calcium isotopes. It 
is the purpose of the present work to close this gap. We 
have chosen in addition to 40 Ca and 48 Ca also to include 
60 Ca, since the latter isotope was studied with a DOM 
extrapolation in Refs. [1, Q. Some preliminary results of 
these FRPA calculations for spectroscopic factors were 
reported in Ref. fl4j but the emphasis in the present work 
is on the properties of the microscopically calculated self- 
energies. The resulting analysis is intended to provide 
a microscopic underpinning of the qualitative features of 
empirical optical potentials. Additional information con- 
cerning the degree and form of the non-locality of both 
the real and imaginary parts of the self-energy will also 
be addressed because it is of importance to assess the 
current local implementations of the DOM method. 

In Sec. Ill Al wc introduce some of the basic properties 
for the analysis of the self-energy. The ingredients of the 
FRPA calculation are presented in Sec. Ill CI The choice 
of model space and realistic nucleon-nucleon (NN) inter- 
action are discussed in Sec. IIIII We present our results 
in Sec. II VI and finally draw conclusions in Sec. [V] 



while at negative energies they give information on fi- 
nal states of the nucleon capture process. Consequently, 
the second term has poles below the Fermi energy (Ep) 
which carry information about the removal of a nucleon 
and therefore clarify the structure of the target state \^q) 
itself. Green's function theory provides a natural frame- 
work for describing physics both above and below the 
Fermi surface in a consistent manner. 

The propagator ([1]) can be obtained as a solution of 
the Dyson equation. 

9<#{E) = 9^(E) + J29 { $(E)Z;s(E)gsp(E), (2) 

in which g(°\E) is the propagator for a free nucleon 
(moving only with its kinetic energy). H*(E) is the irre- 
ducible self-energy and represents the interaction of the 
projectile (ejectile) with the target nucleus. Feshbach, 
developed a formal microscopic theory for the optical po- 
tential already in Ref. [U [22[ by projecting the many- 
body Hamiltonian on the subspace of scattering states. 
It has been proven that if Feshbach's theory is extended 
to a space including states both above and below the 
Fermi surface, the resulting optical potential is exactly 
the irreducible self-energy £*(£) [23| (see also Ref. HI 
and Ref. [25| for a shorter demonstration). 

The above equivalence with the microscopic optical po- 
tential is fundamental for the present study, since the 
available knowledge from calculations based on Green's 
function theory can be used to suggest improvements of 
optical models. In particular, in the DOM, the dispersion 
relation obeyed by is used to reduce the number of 

parameters and to enforce the effects of causality. Thus 
the DOM potentials can also be thought of as a repre- 
sentation of the nucleon self-energy. 



II. FORMALISM 



A. Self-Energy 



In the Lehmann representation, the one-body Green's 
function is given by 
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where a, /3, label a complete orthonormal basis set 
and c a (ct) arc the corresponding second quantization 
destruction (creation) operators. In these definitions, 
l*fe _1 ) arc th G eigenstates, and E* +1 , E^ 1 
the eigenenergies of the (A ± l)-nucleon isotope. The 
structure of Eq. (JTJ) is particularly useful for our pur- 
poses. At positive energies, the residues of the first term, 
(^ / ^ +1 |cJ t |^I / ^), contain the scattering wave functions for 



the clastic collision of a nucleon off the 



ground state, 



For a J = nucleus, all partial waves (£,j,r) are 
decoupled, where £,j label the orbital and total angu- 
lar momentum and r represents its isospin projection. 
The irreducible self-energy in coordinate space (for ei- 
ther a proton or a neutron) can be written in terms of 
the harmonic-oscillator basis used in the FRPA calcula- 
tion, as follows: 



E*(x,x';E) = J2 Zijmtfao-) 



™a ,Tlb 



R n Ar)Kb(E)R nb i{r') 
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where x = r,o~,T. The spin variable is represented by 
ct, n is the principal quantum number of the harmonic 
oscillator, and a = [n a , i, j, r) (note that for a J — nu- 
cleus the self-energy is independent of rrij). The standard 
radial harmonic-oscillator function is denoted by R n e(r), 
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while Xij mj (O, cr) represents the j-couplcd angular-spin 
function. 

We directly calculate the harmonic-oscillator projec- 
tion of the self-energy, which can be written as 



m r a {m r b y 



E — e r ± ir\ 



(4) 



The term with the tilde is the dynamic part of the self- 
energy due to long-range correlations calculated in the 
FRPA, and Y^ h (E) is the correlated Hartrec-Fock term 
which acquires an energy dependence through the en- 
ergy dependence of the G-matrix effective interaction 
(see below). Y^ b {E) is the sum of the strict correlated 
Hartree-Fock diagram (which is energy independent) and 
the dynamical contributions due to short-range interac- 
tions outside the chosen model space. The self-energy can 
be further decomposed in a central (0) and a spin-orbit 
{is) part according to 
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with j>,< = ^±|- The corresponding static terms are de- 
noted by T,^' 1 and and the corresponding dynamic 
terms are denoted by Sq and S| a . 

The FRPA calculation employs a discrete single- 
particle basis in a large model space which results in a 
substantial number of poles in the self-energy Q. Since 
the goal is to compare with optical potentials at positive 
energy, it is appropriate to smooth out these contribu- 
tions by employing a finite width for these poles. We 
note that the optical potential was always intended to 
represent an average smooth behavior of the nucleon self- 
energy @ . In addition, it makes physical sense to at least 
partly represent the escape width of the continuum states 
by this procedure. Finally, further spreading of the in- 
termediate states to more complicated states {3p2h and 
higher excitations that are not included in the present 
calculation) can also be accounted for by this procedure. 
Thus, before comparing to the DOM potentials, the dy- 
namic part of the microscopic self-energy was smoothed 
out using a finite, energy-dependent width for the poles 
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Solving for the real and imaginary parts we obtain 
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where, r implies a sum over both particle and hole states, 
h denotes a sum over the hole states only, and p a sum 
over the particle states only. For the width, the following 
form was used [26f : 



T{E) = - 
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with a=12 McV and &=22.36 MeV. This generates a nar- 
row width near Ep that increases as the energy moves 
away from the Fermi surface, in accordance with obser- 
vations. 

In the DOM representation of the optical potential the 
self-energy is recast in the form of a subtracted dispersion 
relation 
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For the imaginary potential, this is the same as the above 
defined self-energies (01 and it can therefore be directly 
compared to the DOM potential. For the real parts we 
will employ either the normal or the subtracted form in 
the following as appropriate. 



B. Volume Integrals 

In fitting optical potentials, it is usually found that vol- 
ume integrals are better constrained by the experimental 
data JEHU- For this reason, they have been considered as 
a reliable measure of the total strength of a potential. For 
a non-local and ^-dependent potential of the form ^ it is 
convenient to consider separate integrals for each angu- 
lar momentum component, Eg(r, r') and E| s (r, r'). which 
correspond to the square brackets in Eq. (|3]) and decom- 
posed according to ([5]). Labeling the central real part of 
the optical potential with V, and the central imaginary 
part by W, we calculate: 

jf v {E)=4ir [ drr 2 [ dr'r /2 lm £g(r, r'\ E) (11a) 



J^{E)=4n J drr 2 J dr'r' 2 Re S^(r, r'; E). (lib) 

We also employ the volume integral of the central real 
part at the Fermi energy denoted by J F = Jy{E F ), and 
the corresponding averaged quantities 



Jw 9 (E) = J_ J w(E) 

W ££{t} 



Jv 9 (E) = ^- £ J w(E)^ 



(12a) 
(12b) 
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In Eqs. (|12[) . N^y is the number of partial waves included 
in the average and the sum runs over all values of I except 



if otherwise indicated. We also introduce the notation 

The correspondence between the above definitions and 
the volume integrals used for the (local) DOM potential 
in Refs. @, can be obtained by casting a spherical local 
potential U(r) into a non-local form U(r, r') = U(r)S(r — 
r'). Expanding this in spherical harmonics gives 



U(r,r') 



U l (r,r')Yl m (Q.')Y lm (n), (13) 



with the ^-projection 



^( r , r ') = ^(r-rO: 



(14) 



which is actually angular-momentum independent. The 
definition (fTT|) for the volume integrals lead to 

f v = 4tt J dr r 2 J dr'r' 2 U l (r, r') (15) 

= 4tt / U(r)r 2 dr = / U(r) dr , for any I 



and reduces to the usual definition of volume integral for 
local potentials. Thus, Eqs. (fTTj) and (|12 jl can be directly 
compared to the corresponding integrals determined in 
previous studies of the DOM. 



C. Ingredients of the Faddeev- random-phase 
approximation 

The self-energy is shown in terms of Feynman diagrams 
in Fig. [1] The calculations are carried out in two steps 
by following the same procedure as in Ref. [13], where 
further details can be found. First, a configuration space 
is selected that should be as large as possible to account 
for the treatment of nuclear collective motion. We then 
account for the short-range part of a realistic NN interac- 
tion by directly calculating the two-body scattering for 
nucleons that propagate outside the model space. The 
result is the so-called G-matrix that must be employed 
as an energy- dependent effective interaction inside the 
chosen space. The contribution from ladder diagrams 
from outside the model space are also added to the cal- 
culated self-energy and result in an energy-dependent 
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FIG. 1. The self-energy £*(_£?) separates exactly into a 
static (mean-field) term, E°°, and the polarization propaga- 
tors B. (2plh/2hlp) {E) for the 2plh/2hlp motion. These R(E) 
are expanded in terms of particle- vibration couplings as de- 
picted below in Fig. [3] 
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FIG. 2. Expansion of the ph propagator 11(E) in a series 
of ring diagrams. The second line gives examples of time- 
inversion patterns that are generated by the RPA. A similar 
expansion, in terms of ladders diagrams, applies to g XI (E). 
The diagrams are time ordered, with time propagating up- 
ward. 



correction to [see Eq. Q]. When the correspond- 
ing self-energy is calculated, this energy dependence en- 
hances the reduction of the spectroscopic strength of oc- 
cupied orbits by about 10%. A similar depletion is also 
obtained in nuclear-matter calculations with realistic in- 
teractions (Tlj and confirmed by high-energy electron 
scattering data [H, H^. The details^of this partition- 
ing procedure are presented in Ref. 
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discussion, it should be clear that this corresponds to cal- 
culating separately the contribution of propagators that 
lie outside the model space and then to add it to the final 
FRPA results. This does not introduce phenomenologi- 
cal parameters and the calculation should be regarded as 
a microscopic study based only on the original realistic 
interaction. 

In addition to the influence of short-range (and tensor) 
correlations, it is essential to consider the role of long- 
range correlations in which nucleons couple to low-lying 
collective states and giant resonances. This is calculated 
in the second step inside the model space by employing 
the FRPA method. The physics content of the FRPA 
is better summarized by looking at its diagrammatic ex- 
pansion illustrated in Figs. [2] and El The basic ingre- 
dients arc the particle-hole (ph) polarization propagator, 
H a /3,-ys(E), that describes excited states of the A-nucleon 
system, and the two- particle propagator, jS (E), that 
describes the propagation of two added/removed parti- 
cles. These propagators are calculated as summations of 
ring and ladder diagrams in the random-phase approxi- 
mation (RPA). This allows for a proper description of col- 
lective excitations in the giant-resonance region when the 
model space is sufficiently large. The RPA induces time 
orderings as those shown in Fig. [2] for the ph case and 
accounts for the presence of two-particle-two-hole and 
more complicated admixtures in the ground state, which 
are generated by correlations. In FRPA, the R( 2plh \E) 
and R( 2hlp \E) propagators that appear in Fig.[T]are ob- 
tained by rccoupling H(E) and g n (E) to single-particle 
or hole states, as shown in Fig. [3] This is done by solving 
the set of Faddeev equations detailed in Refs. [HI, [l3| . 
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FIG. 3. Left: Example of one of the diagrams for rVv 1k > {E) 
that are summed to all orders by means of the Faddeev 
method. Each of the ellipses represent an infinite sum of 
rings [n(_B)] or ladders [g II {E)\. The diagrams included in 
II(i?) are shown in Fig. [2] and g II (E) is the analogous for 
ladders Right: The corresponding contribution to the 

self-energy obtained from R {2plh) (E) (compare to Fig.©. 



Contributions from ph, particle-particle and hole-hole 
excitations in all possible partial waves are included in 
FRPA as this is required for a complete solution of the 
problem. Moreover, R( 2 P lh ) (E) and R( 2Ii1 p\E) also in- 
clude energy-independent vertex corrections to ensure 
consistency with perturbation theory up to third-order 
to guarantee accurate results at the Fermi surface j3lj . 
We refer the reader to Ref. [TH for more details. 

The reference state employed in calculating the FRPA 
self-energy corresponds to a Slater determinant and is 
chosen to optimally approximate the fully correlated 
propagator ([1} near the Fermi energy. Once the self- 
energy is obtained, a new propagator is calculated by 
solving the Dyson Eq. dp and the full procedure is iter- 
ated to self-consistency [17] . 



III. CALCULATIONS 

Extremely large models spaces are not required for 
the present analysis because we already account for the 
short-range part of the interactions through the parti- 
tioning procedure described in Sec. Ill CI llTI . In the en- 
ergy regime we are interested in, short-range physics af- 
fects mainly the real part of the self-energy in the domain 
of interest. The contributions to the imaginary part are 
not included as they show up at very high positive ener- 
gies which are not considered here |Tl| . The self-energies 
of 40 Ca, 48 Ca and 60 Ca were calculated using the FRPA 
in a harmonic-oscillator model space with frequency Ml 
= 10 MeV. Calculations for 60 Ca were possible in no-core 
model spaces including up to 8 major shells (N max =7) 
and we therefore employed this truncation for all the re- 
sults presented in Sec. IIVI This space is deemed large 
enough to provide a proper description of the physics 
around the Fermi surface and qualitatively good at en- 
ergies in the region of giant-resonance excitations which 
are of interest in this study. 
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FIG. 4. Imaginary volume integral J^y(E) of the 48 Ca self- 
energy calculated with model spaces of different sizes. The 
top (bottom) panels refer to the scattering of a neutron with 
angular momentum l=\ (l=S), Dashed, dot-dashed and full 
lines refer to model spaces of 6, 8, and 10 oscillator shells, 
respectively. These results are for the N3LO interaction. 



Green's function theory — and in particular the 
FRPA — involves infinite summations of linked diagrams. 
This implies that computational requirements scale fa- 
vorably with the increase of the model-space size and 
that the method is size extensive, which allows control- 
ling theoretical errors when increasing the size of the 
system. The FRPA method has been tested in purely 
ab-initio calculations of 4 He in Ref. 141 and was found 



to achieve accuracies comparable to coupled-cluster re- 
sults [32| . The further advantage of the FRPA formalism 
is that it calculates explicitly the effects of all many-body 
excitations including the region of giant resonances. The 
result is a global description of the self-energy over a 
wide range of energies. The FRPA is then the method 
of choice for our purpose of investigating medium-mass 
nuclei in a wide energy domain around the Fermi surface. 
We note that for a calculation of the ground-state energy 
the partition method implies that the contribution of 
high- momentum components still needs to be added [33j j . 
Since these high-momentum components appear outside 
the energy domain of present interest, this issue is of no 
importance here. 

In this work we will focus on averaged properties of the 
self-energy, as described by the volume integrals (fTTI) , for 
which meaningful results can be expected. These will be 
reliable within in a certain energy interval due to the un- 
avoidable truncation of the model space. This window 
is centered around Ep and increases with the size of the 
model space itself. In order to assess these limits, we 
consider the Jw of 48 Ca obtained with the N3LO inter- 
action by Entem and Machleidt [34| for model spaces of 
different sizes. Calculations for this nucleus are possi- 
ble including up to 10 major oscillator shells, as reported 
in 18| . Figure @] shows the proton Jw (E) in the £—0 
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and £=3 partial waves for models spaces of 6, 8, and 
10 shells (and including all orbits with angular momen- 
tum up to £ < 7). As expected, results are similar over a 
range of energies that increases with N max . For higher 
positive energies (in the particle scattering case) Jw is 
expected to increase even further, however, the calculated 
values drop quickly to zero due to the lack of degrees of 
freedom. Based on Fig. @] one can expect that the self- 
energies calculated for N max =7 (8 shells) and discussed 
in this work will be meaningful for energies in the range 
-100 McV< E - E F <100 McV. 

In addition to the chiral N3LO interaction, we have 
performed calculations with the realistic Argonne AV18 
potential [35[ which is representative of a local strongly- 
repulsive core. The J^ 9 for 40 Ca, 48 Ca, and 60 Ca ob- 
tained with the two interactions are compared in Fig. [5] 
The AV18 yields more absorption than the N3LO inter- 
action, for E > Ep, especially in 40 Ca where there is 
about 20% more absorption. Below Ep the absorption 
is only slightly higher. Another important difference is 
that the absorption strength of AV18 is enhanced at en- 
ergies near Ep and on the particle side. Nevertheless, 
it is clear from Fig. [5] that the two interactions generate 
qualitatively similar results. For this reason we will show 
mostly results for the AV18 in the following. 
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FIG. 5. Comparison of the J^ 3 for AV18 (solid line) and 
N3LO (dashed line). 



IV. RESULTS 
A. Angular-Momentum Dependence 

Figures |5]and [7] give an overall example of the features 
of the imaginary (J(y) and real (J v ) part of the self- 
energy. These results are shown for neutrons in 40 Ca, 
employing the AV18 interaction and are separated in par- 
tial waves up to £—5. In Fig. [51 plots of illustrate 
that absorption decreases systematically for increasing £ 
for neutrons, and naturally also for protons. As a conse- 
quence the variation of J v with respect to Jp = J v (Ep) 
obtained from £^°' (Ep), also decreases with increasing 
I (Fig. [7]) on account of the dispersion relation between 
the real and imaginary parts of the self-energy. The effect 
may be partly explained by the truncated model space, 
since the higher ^-channels also have fewer orbits. On 
the other hand, the horizontal lines in Fig. [7l which are 
the contributions of to Jp — J v (Ep), clearly sug- 
gest that most of this decrease must arise from the in- 
dependence implied by the non-locality of the potential. 
Such an ^-dependence suggests that it may be important 
to include non-local features in DOM potentials. For a 
given energy, both the static term and the dynamic term 
have similar radial shapes. In Fig.[5]thc volume integrals 
Jp = J v (Ep) are shown excluding the contribution of 
the dynamic part. Note that because the proton poten- 
tial is not as deep as that of the neutrons, the volume 
integral will be smaller for protons than for neutrons. 
When the calculation is done without the Coulomb po- 
tential, the volume integrals for the protons are compa- 
rable to those for the neutrons. 

This effect of non-locality can be illustrated by taking 
e.g. the energy dependence of the volume contribution 
of a DOM potential [§] and replacing the radial form fac- 
tor by a non-local potential. The radial parameters of 
such a non-local potential employed here correspond to 
the non-local Hartree-Fock potential of Ref. [36| . Such a 
non-local potential is of the form proposed by Perey and 
Buck [37j and contains a gaussian form factor describing 
the non-locality. The results are shown in Fig. Since 
the non-local potential depends on the angle between r 
and r' there is an automatic ^-dependence of the pro- 
jected J^y that exhibit a systematic decrease in absorp- 
tion for increasing £. While it is apparently possible to fit 
clastic scattering data with local potentials, a non-local 
potential has a substantial effect on the interior scatter- 
ing wave function and therefore e.g. on the analysis of 
transfer reactions that rely on such wave functions [38j . 

The possible importance of non-locality for the calcu- 
lation of observables below the Fermi energy was pointed 
out in Ref. [3(|. When the real part of the self-energy 
at the Fermi energy is represented by a truly non-local 
potential, it becomes possible to properly calculate the 
spectral functions below the Fermi energy and observ- 
ables like the charge density. The importance of non- 
locality for the imaginary part of the self-energy sug- 
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FIG. 6. Imaginary volume integral J l w of 40 Ca self-energy for neutrons with I = — 5. 
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FIG. 7. Volume Integrals of Re Eq for neutrons in 40 Ca. The horizontal, dashed lines are the volume integrals of E^°' (Ef) 
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FIG. 8. Angular momentum dependence for the volume In- 
tegrals Jp — Jy(Ef) of E°°' (Ep) excluding the contribution 
of the dynamic part of the self-energy. For each £, results for 
protons are given by solid diamonds and neutrons by solid 
circles. Proton potentials are considerably less attractive due 
to the Coulomb energy. When the Coulomb interaction is 
suppressed (open diamonds) the proton results are close to 
the neutron results. The results shown are for 40 Ca using the 
AV18 interaction. 
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gested by the FRPA calculations may actually provide a 
handle on describing the nuclear charge density for 40 Ca 
more accurately than was possible in Ref. [36j . 

A direct comparison of ^-averaged FRPA volume in- 
tegrals with the corresponding DOM result is made in 
Fig. [TUJ Since the DOM results are calculated from a 
local potential, they must be corrected by the effective 
mass that governs non- locality 0, [3(| , before they can be 
compared with the FRPA results, which are generated 
from non-local potentials. The overall effect of this cor- 
rection is to enhance the absorption. Referring to Fig. 1101 
one can see that the FRPA exhibits different behavior 
above and below Ep than is assumed in the DOM. The 
FRPA predicts that there is significantly less absorption 
below Ep than above, whereas according to the assump- 
tions made in a DOM fit, the absorption is roughly sym- 




E - E P [MeV] 



FIG. 9. Imaginary volume integrals of the volume part of 
a DOM self-energy with a local Woods-Saxon form factor re- 
placed by a non-local form proposed by Perey and Buck. The 
results shown are for £ = (solid), 1=1 (long-dash), 1 = 2 
(long-dot-dash), £ = 3 (short-dash) and I = 4 (short-dot- 
dash). 



FIG. 10. The FRPA results for the average over all ^-channels 
(dashed) are compared with the DOM result (solid) , corrected 
for non-locality. 



metric above and below up to about 50 MeV away from 
Ep 0, ■ While this assumption is made in the local 
version of the DOM, the transition to a non-local imple- 
mentation distorts this assumption of symmetry because 
the attendant correction involving the effective mass is 
different above and below the Fermi energy as can be 
seen in Fig. QUI Since only the absorption above the 
Fermi energy is strongly constrained by clastic scatter- 
ing data, it is encouraging that the ^-averaged FRPA 
result is reasonably close to the DOM fit for both nuclei 
in the domain where the FRPA is expected to be rele- 
vant on account of the size of the chosen model space. 
The simplifying assumptions of a symmetric absorption 
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FIG. 11. Separate partial wave contributions of Jw averaged 
over ^-channels with the same number of harmonic-oscillator 
orbits in the model space. This plot is for neutrons in 40 Ca. 
The dash-double-dotted curve represents the DOM result. 
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TABLE I. Particle-Hole Gaps in MeV 





AV18 


N3L0 


DOM 


Exp. 


40 Ca 


V 


10.7 


12.0 


7.79 


7.23 


7i~ 


7.9 


12.1 


7.20 


7.24 


48 Ca 


V 


4.8 


4.9 


2.83 


4.79 


TV 


11.6 


13.5 


6.78 


6.18 


60 Ca 


V 


4.9 


6.5 


4.95 




TV 


10.4 


12.3 


6.13 





around Ep and locality in the DOM generate unrealis- 
tic occupation of higher ^-values below the Fermi energy 
which is not obtained in the FRPA. More insight into 
this result is obtained in Fig. 1111 where the volume inte- 
gral is averaged over ^-channels with the same number of 
harmonic-oscillator orbits inside the chosen model space 
for neutrons in 40 Ca. Below the Fermi energy the con- 
tribution with 3 and 4 orbits dominate associated with 
the prevalence of low-£ orbits like s,p, and d. Higher i- 
valucs are less relevant below the Fermi energy and this is 
clearly illustrated in the figure. The dash-double-dottcd 
curve illustrates the DOM results also shown in Fig. [TO] 
The DOM result should therefore probably be compared 
below the Fermi energy with curves corresponding to the 
dominant ^-values, whereas above the Fermi energy the 
higher ^-values play a more prominent role. Nevertheless, 
it is clear that the DOM overestimates the absorption 
of partial waves below the Fermi energy that are Pauli 
blocked in agreement with the observations in Ref. |36| . 

Further comparison of FRPA with the DOM self- 
energy is made in Table I for the ph-gap. The AV18 
seems to provide smaller ph-gaps by 1-2 MeV compared 
to N3LO. However, in both cases these gaps substantially 
overestimate the experimental results (see Table I). DOM 
fits from Ref. are also included in the table and are 
typically closer to experiment. 



B. Parity Dependence 

In Fig. Q21 the absorption of the negative parity 
channels is compared with that of the positive par- 
ity channels in 40 Ca, 48 Ca, and 60 Ca. The averages 

(Eevo„ e J w) /^cvon l and (Xodd i J w) l N odd e arc com- 
pared in order to see the trends more clearly. An inter- 
esting feature in 40 Ca is that just below Ep there is more 
negative parity absorption than for even parity, while just 
above Ep the opposite is true. The effect can be under- 
stood in terms of the number of 2plh and 2hlp states, 
which are the configurations beyond the mean-field ap- 
proximation that are closest to Ep. In these states, the 
ph and the hp phonons have negative parity, since the 
holes are in the sd-shell while the particles are in the pf- 
shell. Thus, near Ep, the 2hlp states will have negative 
parity and the 2plh states will have positive parity. 
Proton ph-configurations at low energy continue to 



1 1 i 1 1 1 1 i 1 1 1 1 i 1 M 1 1 i 1 1 1 1 i 1 1 1 1 i 1 1 




FIG. 12. Jw averaged over even ^-channels (solid) is com- 
pared with Jw averaged over odd ^-channels (dashed) 



have negative parity, as the neutron number increases 
in the p/-shell. However, phonons with positive parity 
can be created at energies close to Ep due to the partial 
filling of the neutron p/-shell. So, both parities for 2plh 
and 2hlp states are possible. As a result, in 48 Ca one 
sees little difference between the absorption from nega- 
tive and positive parity states. 

In 60 Ca, which is the next closed shell, the neutron 
p/-shell is filled and the corresponding low-lying neutron 
ph configurations again have negative parity, as in 40 Ca; 
but in this case the neutron holes have negative parity 
corresponding to £ = 1 and 3. Thus, there are more 2hlp 
states with positive parity near Ep for the neutrons. The 
situation for the protons is similar to the case of 40 Ca. 
The inversion of the dominant parity above and below Ep 
is quite general when major shells are filled or depleted 
and also visible in the partial waves separately. 



C. Asymmetry Dependence 

The behavior of the nuclear self-energy with chang- 
ing proton-neutron asymmetry (a = (N — Z)/A) has 
important implications for unstable isotopes. Its un- 
derstanding is fundamental in obtaining proper global 
parametrizations of the DOM so that these can be trusted 
in extrapolations toward the drip lines. Moreover, a 
strong absorption in the optical potential, even if at inter- 
mediate energies, affects the absolute quenching of spec- 
troscopic factors [l8j . Thus, the study of J\y can in prin- 
ciple contribute to the much-debated asymmetry depen- 
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FIG. 13. Asymmetry dependence of the absorption for neutrons and protons and dependence on tensor correlations. The top 
panels shows J^ 9 for 40 Ca (solid) is compared with the results for 48 Ca (dashed), and 60 Ca (dot-dashed). The middle panels 
are obtained by suppressing the tensor component in the AV18 interaction. The difference between the top two panels is plotted 
in the bottom panels to provide a more detailed assessment of the correlations induced by the tensor force. 



dence of spectroscopic factors observed in knockout and 
transfer reactions 0, I38l - l43j . 

The J^ 9 for the three different Ca isotopes are shown 
in the top panels of Fig. [TjJ] These results predict an 
opposite behavior of the protons and neutrons above 
Ep, with the proton (neutron) potential increasing (de- 
creasing) when more neutrons are added, qualitatively 
in agreement with expectations from the Lane potential 
model A recent DOM analysis based on several 

isotopes, including the Ca and Sn chains, employs a sim- 
ilar trend for the volume integrals [l(|. However, the 
same analysis suggests different behavior of the imagi- 
nary surface contributions: neutron surface absorption 
appears to be rather independent of asymmetry, while 
variations are much stronger for protons and for chains 
of isotopes tends to increase with asymmetry. The sep- 
aration between volume and surface effects is an arti- 
fact of the functional form chosen for the optical model 
and such a separation cannot be carried out uniquely in 
a fully microscopic approach like the present FRPA. In 
general, one can argue that most of the physics at scat- 
tering energies below 50 MeV is dominated by surface 
effects which are well-covered by the FRPA, whereas vol- 



ume effects pertain to higher energies, less well-covered 
by the FRPA chosen model space. At energies below 
the Fermi surface, the overall absorption of both proton 
and neutron does not show much variations with chang- 
ing asymmetry. Since the DOM analysis employs less 
data from energies below Ep, this result must be further 
tested in future work. Current DOM implementations 
assume that surface absorption is similar above and be- 
low the Fermi energy, which is clearly not suggested by 
the FRPA results. 

The above pattern, in which one type of nucleon be- 
comes more correlated when increasing the number of its 
isotopic partners, is a rather general feature in nuclear 
systems that is also found for asymmetric nuclconic mat- 
ter (4f| H^. FRPA calculations of stable and drip- line 
nuclei show that this effect results in an asymmetry de- 
pendence of spectroscopic factors similar to that observed 
in knockout reactions, although the overall change from 
drip line to drip line is rather modest [ijj]. We note, 
however, that there also exist other mechanisms that can 
affect this quenching besides the coupling to the giant 
resonance region, including a strong correlation to the ph 
gap [l7j and effects of the continuum at the drip lines (43| . 
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FIG. 14. Effect of the tensor force and charge exchange on 
correlations on the proton- 48 Ca self-energy. The solid curve 
is the imaginary volume integral J^ 3 from the full FRPA 
calculation, while the dashed curve results from removing the 
tensor term in the AV18 interaction. The dash-dot curve is 
obtained by excluding charge exchange from the full calcula- 
tion. The same results are found for neutrons and the other 
Ca isotopes. 



From the characteristics of the above asymmetry de- 
pendence one expects that the nuclear interaction be- 
tween protons and neutrons plays a major role. The ten- 
sor force of the nuclear interaction can provide one such 
mechanism since it is particularly strong in the T = 
sector. Moreover, it has already been shown to influence 
the evolution of single-particle energies at the Fermi sur- 
face [ItJ . To investigate its implication for single-particle 
properties at energies farther removed from Ep, we recal- 
culated J^ 9 by suppressing the tensor component of the 
AV18 interaction. This is shown in Fig.[l4]for protons on 
48 Ca and in the middle panels of Fig. Q2] for all isotopes. 
Its removal results in a drastic reduction of absorption 
at energies \E — Ep\ >30 MeV. Thus, tensor effects give 
an important contribution to scattering at these ener- 
gies. The difference with the complete solution is plotted 
in the bottom panels of Fig. [T3] to highlight the sepa- 
rate effect of the tensor force although this is not unique 
due to the presence of interference terms. It is apparent 
that the tensor force has a very significant effect on the 
correlations far from the fermi surface, but it contributes 
only to the asymmetry dependence of neutron scattering. 
On the other hand, both scattering and negative energy 
states near the Fermi surface are dominated by correla- 
tions other than tensor which thus produce most of the 
asymmetry dependence obtained in the full calculation. 

In Fig. Q3] (dot-dashed line), we have also calculated 
the FRPA self-energy by suppressing charge-exchange ex- 
citations in the polarization propagator H(E). These 
contributions correspond to a mechanism in which the 
proton (neutron) projectile is Pauli exchanged with a 
neutron (proton) in the target and it was argued that 
it could enhance surface absorption due to the presence 
of Gamow- Teller resonances, with strength increasing as 
« 3(N - Z) 0. However, the FRPA results suggest that 



charge- exchange excitations of the target interfere only 
very weakly with the nucleon-nucleus scattering process. 



V. CONCLUSIONS 

In this investigation, an attempt was made to establish 
links between the DOM, an empirical approach to the 
nuclear-many body problem based on the framework of 
the Green's function method and relevant experimental 
data, and the microscopic FRPA approach. An analysis 
of the volume integrals calculated from both approaches 
proved to be a useful link, and on the whole, both the 
DOM and the FRPA produced similar results. However, 
there were some significant and illuminating differences. 

The FRPA exhibits some important shell effects as 
neutrons are added to 40 Ca. In particular, there is a par- 
ity dependence in 40 Ca and 60 Ca, but not in 48 Ca, where 
both parities occur at low energy due to the partial filling 
of the neutron p/-shell. Such an effect has not hitherto 
been taken into account in the DOM. Inspection of the 
imaginary volume integrals generated by the FRPA also 
calls into question the assumption in most DOM analyses 
that the imaginary part is symmetric about Ep for the 
surface absorption. Further insight into the underlying 
physics of DOM potentials is provided by the observation 
that a substantial contribution of the absorption is due 
to the NN tensor force. This influence becomes domi- 
nant at energies about 40 MeV above or below Ep. For 
protons, however, most of the observed asymmetry de- 
pendence of the absorption at positive energy in DOM 
fits appears to be due to central components of the in- 
teraction. For neutrons the decrease in absorption at 
positive energy obtained with the FRPA must be con- 
trasted with the weaker effects deduced so far from DOM 
fits. Noteworthy is also the relevance of the non-locality 
of the absorption process obtained from the microscopic 
FRPA. It leads to an important ^-dependence that may 
play an important role in explaining data like the nuclear 
charge density that are associated with properties of the 
self-energy below the Fermi energy. Its role in scattering 
processes remains to be studied as well and has important 
consequences for the analysis of transfer and knockout 
reactions which are sensitive to interior wave functions 
generated by optical potentials. 
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